The response of laser interferometers to a gravitational wave 
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Laser interferometer detectors are now widely used in an attempt to detect gravitational waves 
(gw). The interaction of the gw with the light circulating in the interferometer is usually explained 
in terms of the motion of the "free" mirrors that form the interferometer arms. It is however 
instructive to show that the same result can be obtained by simply calculating the propagation 
of an electromagnetic plane wave between "free mirrors" in the curved space-time induced by the 
gw. One finds that the plane wave acquires frequency modulation sidebands at the gw frequency, 
as would be expected from the absorption and emission of gravitons from and to the gw. Such 
sidebands are completely equivalent to the time-dependent phase shift imposed on the plane wave, 
that follows from the conventional calculation. 
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I. INTRODUCTION 

In discussing the operation of gravitational wave (gw) 
interferometric detectors [H, 0] it is usually stated that 
the distance between the beam splitter and the mirrors 
at the end of the arms is changed by the gw @, 3 ■ When 
light beams propagate in the arms, the change in the 
proper length of the two arms results in a phase shift 
between the light beams that return to the beam splitter. 
The phase shift is the physical observable that indicates 
the presence of a gw. The above statement is valid if 
the mirrors are free to move (along the axis of the arms) 
and it is expressed in the laboratory frame of reference, 
what is referred to as the Local Lorentz (LL) gauge. It 
is however much easier to calculate the resulting phase 
shift in the Transverse Traceless (TT) gauge. In the TT 
gauge the gw has a very simple form and the coordinates 
of free particles (i.e. of the mirrors) are not changed in 
the presence of the gw, even though their separation docs 
change [!, |(| . 

The geometry of space is defined by the infinitesimal 
interval 01 



ds z 



dx^dx 11 ' 



H, v, — 0, 1, 2, 3, 



(1) 



where summation over repeated indices is implied. For 
weak fields we write the metric tensor as 



gnu = Vfiv + hp V , < 1. (2) 

The flat-space (Minkowski) metric r? Ml/ is chosen to be 
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Consider the interferometer in the x — y plane and that 
a gravitational wave of angular frequency Q is incident 
along the z-axis. In the TT gauge the potential h^ v is 
given by the real part of the following expression 
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h + and h x are the (real) amplitudes of the "parallel" 
and "cross" polarization states of the gw. To be specific 
we will place the origin of the coordinates at the beam 
splitter and orient the interferometer arms along the x- 
and y- axes, as shown in Fig. (1). In the next section 
we calculate, in the TT gauge, the round trip time for 
the propagation of light from the beam splitter to the 
end mirror and back to the origin, in the presence of a 
gw. We will then obtain the same result by considering, 
again in the TT gauge, the direct interaction of the gw 
with the light circulating in the arms. 
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FIG. 1. Configuration of a simple Michelson Interferometer. 
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II. CONVENTIONAL CALCULATION 



We will work in the TT gauge so the coordinates of 
the mirrors are unchanged by the presence of the grav- 
itational wave. We will calculate the time it takes for 
light leaving the beam splitter (the origin) to reach the 
end mirror, where it is reflected, and to return to the the 
beam origin. For light ds 2 = 0, and since the propagation 
is in the x — y plane, Eq. (1) reduces to 



g 00 c 2 dt 2 + gndx 2 + g 2 idy 2 + 2g 12 dxdy = 0. 

For propagation only along the x-axis, and since we can 
set = 0, 
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The time interval for a round trip is obtained by inte- 
grating from the origin to the end of the arm, x = L, 
and back to the origin. If the light leaves at t = to and 
returns at t = t r 



t r — tn 
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At x {t r ) =t r -t = — + frii(t ? .)- Sm f^ T) e iQT . (8) 
c c ill 

For the y-axis we obtain the same result with hu(t) re- 
placed by ti22(t) = — /in(t). Thus the difference in travel 
time between the two arms is 



At(t) = At x {t) - At y [t) = 2[h + e" lUt ] 



, tl Zsin(f2T) einT 



(9) 



A difference in the time of arrival of the light from the 
two arms, implies a phase shift between the two fields. If 
the angular frequency of the carrier is loq , the phase shift 
corresponding to a delay At, is At/> = ujoAt = ckoAt, 
where kg = uj/c is the wavenumber of the carrier. There- 
fore when the gw is normally incident on the interferom- 
eter plane and when the polarization of the gw is along 
the interferometer axes, the observable phase shift is 



iQ(t-T) 



(10) 



where following the convention adopted for Eq. (4), 
Acf>(t) is given by the real part of Eq. (10^. For a more 
detailed discussion of this derivation, see 



where the minus sign before the second integral accounts 
for the fact that on the return trip dt/dx = — \dt/dx\. 
In carrying out the integration we can replace t in the 
exponential by (x/c), because any corrections will be of 
second order in h. We also define the one way travel time 
in the absence of the gw as T — L/c. Thus 
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This suggests that we can write 
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M*r) = h + e- mt " ~ h + e-^ t0+2T \ 

where we have neglected higher order terms in the met- 
ric. As a result, the travel time along the x-axis can be 
expressed as 



III. THE FIELD EQUATIONS IN THE 
PRESENCE OF THE GW 



We will now show that the same result can be obtained 
by considering the direct interaction of the gravitational 
wave with the light propagating in the interferometer 
arms. In this section we find the equations for the elec- 
tromagnetic field in the presence of the gw, and in the 
next section we solve the equations for light propagating 
in the x and/or y arms. We find the presence of (fre- 
quency modulation) sidebands on the carrier displaced 
by the gw frequency f2. This is equivalent to the phase 
shift found in Eq. (10). 

In flat space Maxwell's inhomogencous equations can 
be written in the manifest covariant form as 



d^=f, (ii) 

where the field strength tensor is defined as 

Ffn, = d^A„ - d v A^. (12) 

In a curved space and in the absence of sources, Eq. (JTTJ) 
is replaced by @ 



{V~g g^g vp F Xp ) = 0, 



(13) 



3 



where g = det(g fJiV ) with the field strength tensor still 
defined by Eq. (fT2")l . The dynamical equation (fTB"]) can be 
derived as the Euler-Lagrange equation from the action 
described by the Lagrangian density 



= -\V=9 ff"V p (d»A v - d v A^) (d\A p - d p A x ) . 

(14) 

We will again work in the TT gauge and take the gw to 
be the same as before. In the TT gauge, in the weak-field 
approximation y/—g = 1, and keeping only terms linear 
in hp U , Eq. (fl"3"j) reads 



d^ v - 8^ X F X ") - d^F%) = 0. 



(15) 



The first term is the same as Eq. (fTTj) with j v = 0, 
and describes the propagation of the free electromagnetic 
field. The next two terms act as sources that give rise to 
new electromagnetic fields generated by the interaction 
of the gw with the field F^ v of the light circulating in 
the arms. For the choice of gw (normal incidence toward 
the negative z-axis), expansion of Eq. (14) in E and B 
field components yields the four equations 

v = : V • E - h+e-^ m+k ^ (d x E x - d y E y ) 
- h x e-^ nt+k ° z \d x E y + d v E x ) = 0, 

v = 1 : d ct E x - (V x B) x + h+e~' l{nt+k ^ 

x (iQE x + iknBy - d ct E x - d z B y ) 
+ h x e- lint+k " z} (iQEy + ik n B x - d ct Ey + d z B x ) 
= 0, 

v = 2 : d ct Ey - (V x B)„ - h + e-' l(nt+kaz) 

x (iflEy + ik n B x + d z B x - d ct E y ) 
- h x e- l(nt+knz) {inE x + ik n By - d ct E x - 8 z B y ) 
= 0, 

v = 3 : d ct E z - (V x B), + h+e^ nt+k ^ (d x B y + d y B x ) 

(d X B X - dyBy) = 0. (16) 



IV. SOLUTION OF THE EQUATIONS OF 
MOTION 



The electromagnetic (em) fields E, B in the arms con- 
sist of the carrier field Eq, Bo due to the external source 
(the laser beam), and the additional field E + ,B + , gen- 
erated by the source terms in Eq. (|T7|) . 
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Since Q <C ujq we can drop the terms in O and fcn. 
We also make use of the condition z = and choose 
h-\- 7^ 0, hx = 0. This leads to the simplified equations 
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(17) 

In the next section we will show how to solve these equa- 
tions when a plane (light) wave of angular frequency ujq 
propagates along the x and/or y axis. 
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B = B, 
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(18) 



We describe the carrier as a plane wave propagating along 
the a;-axis and polarized along the z-axis, and of constant 
amplitude Aq 
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To the sideband fields E+,B+, we assign slowly time- 
dependent amplitudes E+(t), B + (t) and a phase factor 
{u)+t — kox) where uj + = ujq + O, namely, 



E^ 



E+{t)e 
B+(t)e 
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(20) 



We will solve the equations in (fT"7|) perturbatively not- 
ing that E + ,B + are of order h with respect to E ,B - 
Thus on the left side of Eq. (|17|) we must use the full 
fields E, B but on the right side it suffices to retain only 
E ,B . With the choice of Eq. (fl9|) for the carrier field 
only the u = 3 equation is nontrivial; the first three equa- 
tions are satisfied automatically with the choices made 
in Eqs. (19,20). Writing out the part of order h, for the 
v = 3 equation in (17), gives 
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The remaining terms of the equation 
dE 0z dB^y 
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(21) 



(22) 



describe the propagation of the free carrier field. The 
fields, both carrier and sideband satisfy Maxwell's dual 
equation, namely d^F^ = 0. For our choice of direction 
of propagation and polarization, this condition reduces 
to 



dE 



+z 



dB 



+u 



0. 



(23) 



Combining Eqs. 
for E +z ,B +y 



dx cdt 
2Tj) and (|23|) leads to wave equations 
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= -h + k + k A e- l( - u '+ t - kox \ 

+h+k 2 A e- l( - UJ + t - k ° x \ (24) 



where we have introduced k + = lo + / c. 

We now use Eqs. (|2"0)) and keep in Eqs. (24) only the 
first time derivatives dE + (t) / dt, dB + (t) / dt of the ampli- 
tudes, since E + (t), B + (t) vary slowly. We can also ap- 
proximate oj\jc 2 — fcg ~ 2fiwo/c 2 and find the following 
two eqs for the slowly varying amplitudes 



dE^ 
~dt 
dB + 
dt 



- i — IIE+ = i . 

cj + 2 



Wn„„ 

i— VlB^ 



(25) 



with solutions 



A = An 



iA 4> 



A+ 



A e 



i<j>A 4 



(31) 



which shows that the phase of the carrier is shifted by 
exactly the same amount as was calculated in section 2. 
When the carrier propagates in the y-arm the sign of h 
and thus also of A + (2T) is reversed; after subtracting 
the fields returning from the two arms, we find for their 
relative phase difference 24> A +, as we had obtained in 
Eq. (10). 

Consider now a gw that depends sinusoidally on time, 
which wc express at z = by 



h+(t) — h + cos fit = h + 



(32) 



It is clear that in this case both an upper and lower side- 
band will be present. The fields are 



(26) 



where we assumed the initial conditions E + {t = 0) = 
B + (t = 0) = 0. Therefore, to the approximation 57 <C w 
that we are using, the sideband fields are 



E, 



iE, 



OVA- i( w _t-fe x) 



where 



E_ = 

(33) 

2T ; to see this change 



OA- = 4> A + = PAe 
fi ->• - ft in Eq.(30). At t = 0, x = the two sideband 
fields have equal imaginary parts and opposite real parts. 
This is the condition that corresponds to phase (or fre- 
quency) modulation when the sidebands are combined 
with the carrier 
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(27) 



The amplitudes E+(t), B + (t) of the sideband fields are 
zero at t = and grow in time as a parametric amplifi- 
cation process. We will designate them by A + (t). The 
maximum value of A+(t) is determined by the losses in 
the arm cavities, however the carrier amplitude, Aq, re- 
mains constant. For one round trip t = 2L/c = 2T and 
the sideband amplitude has the value 



A + (2T) = ih + Lk A e 



iQT 



sin nr 
SIT 



(28) 



We see that after one round trip the amplitude Aq of 
the carrier acquires a small complex part 



A+(2T)=iA (t> A+ , 



(29) 



where 



b A e m = h+ k L S -^fe^ T « 1. (30) 



Thus the amplitude of the propagating light can be writ- 
ten as 



E = Ene~ i ^° t ~ koX ^ 
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(34) 



In conclusion a sinusoidal gw interacting with a plane 
wave carrier (under the appropriate geometry) imposes 
upper and lower sidebands, or equivalcntly contributes a 
time-dependent phase shift. We obtained this result by 
considering the propagation of the carrier in the space- 
time of the gw, without referring to the change in the 
round-trip time of travel to and from the end mirror. 



V. DISCUSSION 

Our derivation is based on a subtle point, the use of 
the TT gauge. We can use the TT gauge only because the 
end points of the light travel (the mirrors) are free. If the 
mirrors were fixed we would have to do the calculation in 
the LL gauge, i.e. in the laboratory frame of reference. In 
the LL gauge fl(il - [T3l ] the amplitude of the gw is modified 
from the form given by Eq. (4), the relevant part of the 
metric being 
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h m = —h(t)(x 2 -y 2 ). (35) 

For low frequency gw's the resulting phase shift 
when the mirrors are fixed is much smaller than for 
free mirrors, and this is why interferometric gw de- 
tectors are constructed with suspended mirrors. For 
the details of the calculation in the LL gauge see ref. [HI ■ 

The presence of sidebands has a direct physical 
interpretation in terms of absorption and stimulated 
emission of gravitons from/to the gw. Since the gw field 
is highly classical, [the occupation number for h ~ 10~ 23 
and / - 100 Hz, is n = iV g /(A/27r) 3 - 10 33 ], both 
processes have the same probability. Absorption leads 
to the upper sideband, while emission to the lower one. 
There is no energy exchange between the optical and 
gravitational fields, but only a phase shift. 

In our analysis we have assumed that the carrier is a 
plane wave propagating toward the positive x-direction, 
see Eqs. (19). In practice the carrier is reflected by a 
mirror at x = a, thus the Electric field should vanish at 
x = a. This condition is satisfied by writing the fields in 
the form 



If a mirror is also placed at x = to reflect the carrier 
and koa = mr, a standing wave is established in the 
region 5= x ^ a. This does not modify the conclusions 
that we have reached. The case when the carrier is a 
standing wave is treated by Cooperstock and Faraoni 
(13 ] and leads to the same results as in Eq. (34). The 
interaction of an electromagnetic and gravitational field 
was first discussed by Gertsenshtein and Pustovoit [l5[ 
in 1963, and subsequently by others fl2l [l6l - [l8| . 

We have also shown that a gw will couple to the car- 
rier in a single arm. The interferometer configuration has 
been chosen because it is technically advantageous. The 
light returning from the two arms is adjusted to interfere 
destructively at the detection point in the absence of a 
gw. Thus when the gw induces a phase shift, a signal ap- 
pears over a null background (excluding noise). By using 
multiple traversals in the arms one can increase the effec- 
tive length, and thus the phase shift, significantly. The 
main limitation in effective arm length (apart from losses 
in the optics) is related to the frequency of the gw. We 
have seen in Eq. (10) that the phase shift is modulated 
by the "form factor" sin(QT)/QT, and therefore fLL/c 
must be kept small. 
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